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ABSTRACT 

We have analysed a sample of satellite and primary galaxies in the 2dF galaxy redshift 
survey. In our study we find a strong statistical evidence of the Holmberg effect (that 
is a tendency for satellites to avoid regions along the line defined by the primary 
plane) within 500 kpc of projected distance to the primary. This effect is present only 
when we restrict to objects with radial velocity relative to the primary \Av\ < 160 
km/s which correspond approximately to the mean of the distribution. We explore the 
dependence of this anisotropy on spectral type (77) , colours, and luminosities of both 
primaries and satellites, finding that objects with a low present-day star formation 
rate present the most significant effect. 
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1 INTRODUCTION 



The origin and evolution of structures in the universe can 
be suitably understood within the hierarchical scenario. In 
such schemes, present-day systems are the result of the ag- 
gregation of several sub-structures, and in particular, satel- 
lite galaxies can contribute significantly as building blocks of 
galaxy formation. However, present-day systems of satellites 
have evolved since their formation; and therefore the obser- 
vations of such systems can provide invaluable information 
on the formation and evolution of galaxies. 

There are several works that explore the characteristics 
of systems of satellites. Holmberg (1969) found that the dis- 
tribution of galaxy satellites presents a lack of objects within 
30° of the primary disk plane (0 < 30°). He searched for 
companions of nearby bright disk galaxies (the primaries) 
within 50 kpc of projected separation (rp). His sample was 
composed of 218 satellites including both physical systems 
and interlopers that contaminated his photometrically se- 
lected sample. He found 45 objects within <j> < 30° while 
173 had (p > 30°. If the position angles relative to the pri- 
mary plane are homogeneously distributed, then the number 
expected within (p < 30° would be ~ 73. After removing con- 
tamination by presumed interlopers (which are expected to 
be distributed isotropically) he concluded that there were 
no physical objects with (j> < 30° . 

The results remained controversial since Valtonen et al. 
(1978) found an excess of compact satellites with small an- 
gles relatives to the primary disk. Based on this work, Byrd 
et al. (1987) argued that the anisotropy of the Holmberg 



sample was mainly due to selection effects for he may have 
not considered these compact dwarf objects. The authors 
proposed a dynamical stripping of the gas content of the 
satellite as a possible explanation for the excess of such ob- 
jects near the primary plane. 

Some years later, the problem was investigated again 
in a statistical way by Zaritsky et al. (1997) using redshift 
data. They selected a sample of 115 satellites of 69 nearby 
primaries, these satellites have a maximum projected sepa- 
ration Vp < 500 kpc. These authors could not confirm the 
Holmberg effect in the inner 50 kpc since they had few ob- 
jects at this distance , but they found a similar anisotropy 
effect with reliable statistical significance at larger separa- 
tions, Vp > 250 kpc. 

One simple starting point to approach this problem is to 
assume an isotropic initial distribution of satellites and that 
evolution generates the observed anisotropics. This scenario 
has been explored by Quinn & Goodman (1986), who found 
that anisotropy in the distribution of satellite positions is not 
efficiently produced by dynamical evolution, even at small 
separation, Vp < 50 kpc. 

On the other hand, numerical simulations like those car- 
ried out by Abadi et al. (2002) indicate that, in a hierar- 
chical scenario, the formation of disks involve substantial 
accretion of satellites which strongly infiuence the angular 
momentum vector and therefore the primary plane. Thus, 
observed anisotropies in the distribution of the remaining 
satellites could be a direct consequence of this process. In 
addition to this work, Pefiarrubia et al. (2002) also based 
on numerical simulations, established that the time interval 
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before the satellite disruption strongly depends on the dark 
matter halo shape, the initial orbital eccentricities and in- 
clinations, and on the mass fraction of satellites relative to 
the primary. 

Our work aims to search for anisotropies of the angu- 
lar distribution of satellites relative to the primary plane as 
well as exploring the possible dependence on satellites prop- 
erties. Our analysis is based on the study of an unprecedent- 
edly large statistical sample extracted from the 2dF galaxy 
redshift survey. In section 2 we present the data and selec- 
tion criteria, in section 3 we show the main results, and in 
section 4 we summarise our conclusions. 



2 PRIMARY AND SATELLITE SAMPLES 

The largest complete redshift survey available to date is the 
2dFGRS which comprises information of objects on approx- 
imately 1500 square degrees of the sky in the north (NGC) 
and south (SGC) strips, and in several random fields. This 
catalogue has been obtained with the multifibre technology 
at the Anglo- Australian Observatory, and provides photo- 
metric and spectroscopic data for 232155 galaxies with mag- 
nitudes bj < 19.45. It contains positions, magnitudes (in bj 
and R band), eccentricities, position angles, redshifts (well 
determined for 221414 galaxies) and the spectral classifica- 
tion through the rj parameter (which has a high confidence 
for objects with z < 0.1). 

The selection criteria adopted in this work are the usual 
in the literature. From 2dFGRS we have selected 1567 pri- 
mary galaxies with z < 0.1 restricted to absolute magnitude 
Bj < —18. We have imposed an isolation criterion by consid- 
ering primaries with no companion galaxies within projected 
distance < 700 kpc and relative radial velocity difference 
|Ai;| < 1000 km/s, and with absolute magnitude difference 
j^comp _ j^prtrn ^ ^ The value of the Hubble constant in 

this work is assumed to be Ho = lOO/i km/s Mpc~^ . 

After the selection of the primary sample, we searched 
for satellite galaxies around each primary. We define as satel- 
lites galaxies with rp < 500 kpc and |A«| < 500 km/s 
which are at least 2 magnitudes fainter than their primary, 
j^prim _ j^sat y 2 With all these restriction our final sample 
consists of 3079 satellites. 

A few primaries have too many satellites (up to A" = 19 
objects). These objects are likely to be contaminated by 
group-like systems with the primary as the brightest group 
member. In order to exclude these spurious systems we have 
restricted the sample to 1498 primaries with N < 5 out 
of the 1567 original primaries. Furthermore, we notice that 
the removed systems show systematically high mean velocity 
dispersion, close to typical value of groups of galaxies, sup- 
porting our conservative cut-off in the number of members 
per system. 

We have also imposed an additional restriction to the 
eccentricities of the primaries in order to deal with well de- 
termined relative angles <f) of satellite position with respect 
to the primary plane projected into the first quadrant. Thus, 
we have not considered primaries with e < 0.1 
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Figure 1. Distribution of position angles of satellites for two 
subsamples corresponding to high and low relative velocity to the 
primary. The solid line represents satellites with | Ai)| > 160 km/s 
and the dotted line satellites with |Ad| < 160 km/s. The adopted 
threshold \Av\ = 160 km/s is approximately the mean value of 
the satellite relative velocity distribution (163.9 km/s). 



3 RESULTS 

The sample selected under the criteria described in the pre- 
vious section allows us to investigate the satellite position 
angle distribution relative to the primary plane disk (0). 

We have divided the sample into objects with high or 
low \Av\ compared to \Av\ = 160km/s (a value close to 
the mean |Aw| = 163.9 km/s). The distribution of satellite 
position angles for these two subsamples are shown in Fig. 
Q From this figure one can notice the deficiency of satellites 
with low (f) values in the subsample of satellites with relative 
velocity |A| < 160 km/s. On the contrary, the high velocity 
subsample is roughly isotropic. 

We may quantify the departure from isotropy defining 
the ratio / of objects in planar {(f) < 30°) to polar (0 > 60°) 
angular positions / — A''<3o/A>6o, for isotropy / = 1. The 
resulting values are: / = 0.80 H — 0.04 for \Av\ < 160 km/s 
and / = 0.98 -I- -0.08 for \Av\ > 160 km/s. Quoted uncer- 
tainties in this work were calculated based on 20 bootstrap 
re-samplings of the data. 

Other useful statistical measures of this anisotropy ef- 
fect can be obtained by fitting a function to the (j) distribu- 
tion. We have adopted a double cosine function N/N — 
Acos{2(l}) + B, and a linear fit N/N — aunfj) + b where 
iV is the mean number of objects and A and aun are the 
anisotropy parameter in each fit. For an isotropic distribu- 
tion on (f) we expect A — aun ~ 0. Either function, the 
double cosine and the linear relation, are suitable to fit the 
data, see for instance Fig. |21 We notice, however that the 
double cosine function is less sensitive to noise fluctuations 
at the extremes due to poor number statistics and therefore 
might provide more conservative estimate of the deviation 
from isotropy. For this reason we will mainly refer to the A 
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Table 1. Anisotropy parameters aii„ and A, and their bootstrap 
derived error estimate. 



Sample 


Characteristics 


N 






A 




S — vel^ 


\Av\ < 160 km/s 


0.4 


0.1 


1276 


0.12 


0.04 


S — vel^ 


\Av\ > 160 km/s 


0.0 


0.1 


916 


0.01 


0.04 


S - rp< 


\Av\ < 160 km/s, 














rp < 250 kpc h-^ 


0.2 


0.1 


713 


0.11 


0.05 


S - rpy 


\Av\ < 160 km/s, 














rp > 250 kpc h-^ 


0.2 


0.1 


691 


0.12 


0.06 



parameter although both A and ann values are listed in the 
tables corresponding to the different subsamples analysed. 

In table we list the number of satellites in each sub- 
sample (|At;| < 160 km/s and \Av\ > 160 km/s) and the 
corresponding aun and A values. The A parameter for sub- 
sample S — iie/< departs from isotropy by more than ^ 3a 
which contrasts with the behaviour of subsample S — iieZ> 
where A ~ 0. This reflects the trends shown in Fig. Q and 
is consistent with our previous estimate of /. The fact that 
the anisotropy signal is not significant in the high relative 
velocity subsample {\Av\ > 160 km/s) suggests a higher 
contamination due to interlopers for these objects (see for 
instance Zaritsky et al. 1993). 

The number of objects per bin is A'^ = Ns{4>) + Ni{(j)), 
where Ns and Ni are the number of physical satellites and 
interlopers as function of angle, respectively. False satel- 
lites should have a homogeneous distribution of (p angles, 
so we take Ni{(j>) = A'^;. On the other hand, the distri- 
bution of real satellites respect to the primary plane may 
present anisotropies. The observed anisotropy parameter 
Aobs is related to the real anisotropy parameter Areai by: 
Aobs = {< Ns > / < Ns > +Ni)Areai- Thus, cveu a modest 
observed anisotropy ~ 10%, can imply a substantially large 
real anisotropy. 

As suggested by Zaritsky (1992), the relative fraction 
P/T of objects with At; > for a given sample can provide 
a useful measure of the degree of contamination given the 
different volumes involved. If a sample is seriously contami- 
nated by interlopers, we should obtain P/T » 0.5. 

Therefore, we calculated P/T for two subsamples, satel- 
lites with |Au| < 160 km/s and satellites with \Av\ > 160 
km/s. For these two subsamples P/T values are 0.54 ± 0.02 
and 0.51 ±0.03 respectively, so that from this particular test 
we cannot argue for substantial contamination in the high 
velocity subsample. Thus the different behaviour of high and 
low relative velocity satellites could be a real effect. 

In his pioneer work, Holmberg found an excess of po- 
lar satellites within 50 kpc from the primary; neverthe- 
less, Zaritsky et al. (1996) could not reproduce this result, 
but they found such an excess at larger separations (250 
kpc< rp <500 kpc) from the primaries. We have searched 
for possible dependences of the detected anisotropies on pro- 
jected separations to the primaries. The resulting anisotropy 
parameter A are displayed in table where it can be ap- 
preciated the lack of strong dependence on projected sepa- 
ration. From now on, we restrict our sample of satellites to 
those objects with |Ai;| < 160 km/s. 

In order to investigate if there are any preferred satel- 
lite characteristics which show a more significant anisotropy 
signal, we have analysed the distributions of absolute mag- 
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Figure 2. Distribution of satellite position angles for subsample 
with |Aj;j < 160 km/s. We show two possible fit to the data. In 
solid curve is the linear relation and in long-dashed ones is the 
double cosine function. 



Table 2. The anisotropy parameters 
different satellites properties. 



and A for samples of 



Sample 
SI 
S2 

S3 
S4 

S5 
S6 



Characteristics 

gsat ^ 3 



< 1.1 
> 1.1 



< 0.83 
> 0.83 



N 






A 




633 


0.3 


0.1 


0.15 


0.04 


643 


0.1 


0.1 


0.09 


0.04 


580 


0.3 


0.1 


0.16 


0.04 


643 


0.1 


0.1 


0.08 


0.04 


644 


0.3 


0.09 


0.15 


0.05 


632 


0.2 


0.09 


0.09 


0.06 



nitude and spectral type of satellites in polar (0 > 60°) and 
planar (<^ < 30°) positions. The objects with tj) < 30° and 
(j} > 60° are 385 and 477 respectively. We show in Figl^a) 
and b), the corresponding distributions of absolute magni- 
tude Bj and spectral type parameter rj for these two cases. 
The error bars correspond to Poisson uncertainties. It can be 
appreciated in these figures that the excess of objects with 
4> > 60°, is distributed quite uniformly along the explored 
range of luminosities and spectral types. 

In order to quantify the anisotropy for these subsamples 
we show in table |5| the corresponding values of aun and A 
parameters. 

The adopted thresholds in absolute magnitude Bj , spec- 
tral type ri and colour index bj — R were chosen in order to 
divide the samples into approximately equal number sub- 
sets. 

By inspection to table 2 it can be seen that the 
anisotropy of sample S2, although less statistically signifi- 
cant that in sample SI, still shows an excess of polar satel- 
lites at a ^ 2.5a level. Thus, this result does not support 
the hypothesis of Byrd et al. (1987) where the lack of satel- 
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Figure 3. Different satellite parameter distributions. Dashed 
(solid) curves are polar (coplanar) angular positions. In upper 
panel is Bj distribution, in the bottom is represented the spec- 
tral type r\ of the satellites. The error bars correspond to Poisson 
uncertainties. 
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Figure 4. Double cosine fits for different subsamplcs: a) 5 — fe/< , 
h) S3, c) P2 and d) P3. 



Table 3. The anisotropy parameters a;i„ and A for samples of 
different primary properties. 



Sample 


Characteristics 


N 




O-a 


A 




PI 




630 


0.2 


0.1 


0.10 


0.07 


P2 




643 


0.3 


0.1 


0.14 


0.05 


PS 


if < -1.4 


686 


0.4 


0.1 


0.18 


0.03 


P4 


•qP > -1.4 


581 


0.1 


0.1 


0.03 


0.06 


P5 


{bj - R)P < 1.13 


645 


0.1 


0.1 


0.08 


0.05 


P6 


{bj - R)P > 1.13 


631 


0.3 


0.1 


0.16 


0.05 



lites with (j> < 30° in Holmberg's data would correspond 
to a combination of dynamical stripping of gas and selec- 
tion effects, where faint compact satellites should be found 
preferentially near the primary plane. 

The results for subsample S3 indicate that the satel- 
lites with low star formation activity present the strongest 
anisotropy signal compared to those of high star formation 
activity (subsample S4). Subsamples S5 and S6 indicate that 
the satellites characterized by low colour index values show 
a stronger excess near the primary poles. This can be appre- 
ciated from Fig 13)), where is evident that the intermediate 
spectral type satellites ( — 1.4 < rj < 1.1) show a strong defi- 
ciency at < 30°. 

We have also explored for possible dependencies of the 
anisotropy on properties of the primaries, namely luminos- 
ity, spectral type and colour index. In tablej^we list N, aun, 
a a, A and a a for subsamples selected according to primary 
properties (regardless of satellites characteristics). 

Different ranges of absolute magnitudes of the primaries 
show different (j) distributions (see first two lines in table 



13 . For faint primaries (P2) the angular positions of satel- 
lites are well reproduced by our double cosine fit, where the 
anisotropy parameter differs from zero by ~ 3a. However, 
for bright primaries (PI) this effect is not statistical signifi- 
cant. In a recent paper, Prada et al. (2003) have shown that 
the luminosity and halo mass of a galaxy are closely related. 
This results suggest that the anisotropy of the satellite dis- 
tribution would not have a tidal origin, since in this case 
we would expect satellites of bright primaries to show the 
strongest anisotropy, contrary to our findings. 

On the other hand, we detect a tight dependence of the 
satellite anisotropy distribution on the spectral type of the 
primary. In table |3] it can be seen that satellites of passive 
star forming primaries (P3) show a departure from isotropy 
of at least ~ 6(7. A different behavior is found for the sam- 
ple of primaries with significant star formation activity (P4) 
which show a remarkable isotropic distribution of satellites. 
Regarding the colour indexes of primaries we find a more 
significant effect for primaries with larger values of bj — R 
colour index (sample P6). 

In Fig. |1| we show four different subsamples with the 
high anisotropy signal; a) satellites with |Au| < 160 km/s 
(5* — vel^), b) satellites with rj" < 1.1 (S3), c) and d) 
primaries with absolute magnitude Bj > —20 (P2) and 
rjP < —1.4 (P3) respectively. 

A useful test to check the real presence of the signal 
in our analysis is to restrict to primaries with high eccen- 
tricity values (e). We have performed this test for subsam- 
ple P3 (primaries with low star formation activity) which 
show a large anisotropy effect. In table |1] we show the re- 
sults obtained for this sample for different threshold values 
of e. It can be seen that the strongest signal corresponds 
to primaries which are nearly edge-on (e > 0.4) which gives 
substantial confidence to our findings. 
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Table 4. The anisotropy parameters aii„ and A for samples of 
passively star forming primaries (rj^ < —1.4) with different ec- 
centricities. 



Sample 


Characteristics 


N 




0"a 


A 


aA2 


P3 


e > 0.1 


686 


0.4 


0.1 


0.18 


0.03 




0.1 < e < 0.4 


338 


0.2 


0.1 


0.15 


0.08 


P3(, 


e > 0.4 


337 


0.3 


0.1 


0.26 


0.08 



4 SUMMARY 

We have carried out an analysis of the anisotropy of satel- 
lite angular positions with respect to the primary plane. 
The sample was selected in redshift space from the 2dF- 
GRS and is large enough for studying possible dependencies 
of anisotropies on both primary and satellite properties. 

We confirm the existence of the effect first claimed by 
Holmberg, that is, the preference of the satellites to lie near 
the pole of the primary. We note however that this effect is 
present up to relative projected separations rp < 500 kpc. 

We use an anisotropy parameter A suitable to measure 
the departure from isotropy in the different subsamples anal- 
ysed. Nevertheless, this parameter indicates a lower limit to 
the real anisotropy, given the contamination by false satel- 
lites expected to be uniformly distributed with respect to 
the primary plane. 

The deviation from isotropy is detected at least at 
~ 3cr level for objects with low values of relative velocity 
Aul < 160 km/s while the high velocity subsample is con- 
sistent with isotropy. Although in principle this effect could 
be related to different levels of contamination by interlopers, 
Zaritsky 1992 test provides comparable values of contami- 
nation in both high and low relative velocity subsamples. 

By inspection to the luminosity of the primaries, we 
conclude that the anisotropy of the distribution of satellites 
around faint primaries shows a high statistical significance 
(^3a). 

The most important signal was detected in the sample 
of satellites corresponding to passive star forming primaries, 
where the distribution of (j) angles departs from isotropy at 
the ~ 6o" level; while the opposite behaviour is present in 
primaries with rj > —1.4, for which the A parameter results 
consistent with zero. Moreover, we also notice that the signal 
is even higher when we raise the eccentricity threshold of 
the primaries. We also have considered the different spectral 
types of the satellites, finding the strongest anisotropy for 
satellites with 77 < 1.1 (objects with little present day star 
formation activity). We notice that primaries with a low 
value of r] in our sample also have satellites with poor star 
formation activity. 

The isolation criteria imposed to our sample lowers the 
probability of including a significant number of elliptical 
galaxies. Thus, most primaries in our sample are expected 
to be of late type morphology. Due to the morphology den- 
sity relation we expect only 10% of our isolated sample of 
primaries to be ellipticals (Whitmore et al. 1993). 

We could in principle interpret the Holmberg effect 
as being mainly caused by evolutionary processes. Systems 
dominated by old stellar population objects have formed 
early, so the physical processes that may have removed satel- 
lites near the primary plane have been present for a large 
time. By contrast, in younger primary and satellite systems, 
we would expect a larger fraction of objects that have not 



been accreted yet on to the primary so that they would show 
an isotropic distribution of satellites. However, as discussed 
by Quinn & Goodman (1986), it may be difficult to pro- 
duce anisotropies in the distribution of satellites from pure 
dynamical evolution, specially at large scale rp < 500 kpc. 
More hkely, we suggest that according to Abadi et al. (2002) 
results, such anisotropies could be the remnants of the for- 
mation of the disks of the primaries, where substantial or- 
bital momentum is transfered into primary rotation during 
the early stages of galaxy formation. Thus, the anisotropy 
effects we detect at large projected separations, could be the 
relics of such mechanism. 

We also argue that the significant dependence of the 
anisotropy effect on the relative velocity between the satel- 
lite and the primary could be an important point of dis- 
agreement between the results of previous works. 
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